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ABSTRACT 


Recently  there  has  been  much  work  devoted  to  considering  some  of  the  many  and  varied 
interaction  mechanisms  which  may  be  operative  in  three-dimensional  boundary  layer  flows. 
Here  we  are  concerned  with  resonant  triads  of  crossflow  vortices.  In  contrast  to  much  of  the 
previous  work  we  examine  the  effects  of  interactions  upon  resonant  triads  where  each  member 
of  the  triad  has  the  property  of  being  linearly  neutrally  stable;  then  the  importance  of  the 
interplay  between  modes  can  be  relatively  easily  assessed.  We  concentrate  on  investigating 
modes  within  the  boundary  layer  flow  above  a  rotating  disc;  this  choice  is  motivated  by 
the  similarity  between  this  disc  flow  and  many  important  practical  flows  and,  secondly,  our 
selected  flow  is  an  exact  solution  of  the  Navier-Stokes  equations  which  makes  its  theoretical 
analysis  especially  attractive.  Firstly  we  demonstrate  that  the  desired  triads  of  linearly 
neutrally  stable  modes  can  exist  within  the  chosen  boundary  layer  flow  and  then  subsequently 
obtain  evolution  equations  to  describe  the  development  of  the  amplitudes  of  these  modes  once 
the  interaction  mechanism  is  accounted  for.  It  is  found  that  the  coefficients  of  the  interaction 
terms  within  the  evolution  equations  are,  in  general,  given  by  quite  intricate  expressions 
although  some  elementary  numerical  work  shows  that  the  evaluation  of  these  coefficients 
is  practicable.  The  basis  of  our  work  lends  itself  to  generalisation  to  more  complicated 
boundary  layers  and  effects  of  detuning  or  non-parallelism  could  be  provided  for  within  the 
asymptotic  framework. 

^Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract 
No.  NASl-18605  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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§1.  Introduction. 


The  desire  to  understand  the  phenomenon  of  transition  to  turbulence  in  fluid  flows 
has  motivated  a  wealth  of  studies  concerned  with  the  instability  of  three-dimensional 
boundary  layer  flows.  Here  we  address  the  problem  of  providing  a  self-consistent 
asymptotic  description  of  the  interactions  within  a  resonant  triad  of  travelling  cross- 
flow  vortices  occurring  in  a  three-dimensional  boundary  layer  at  a  large  Reynolds 
number.  In  particular,  we  relate  our  analysis  to  vortices  in  the  boundary  layer  of  the 
flow  induced  by  a  rotating  disc.  This  flow  is  susceptible  to  instabilities  similar  to  those 
which  occur  in  flows  over  a  swept  wing;  a  situation  which  has  practical  relevance  to 
the  development  of  Laminar  Flow  Control  wings.  Further,  there  is  an  exact  solution 
of  the  Navier  Stokes  equations  which  describes  the  rotating  disc  flow  which  makes 
this  particular  flow  especially  suitable  for  a  theoretical  analysis  of  the  interaction  of 
travelling  cross-flow  vortices. 

The  resonant  interaction  of  a  triad  of  waves  is  one  of  a  number  of  possible  mecha¬ 
nisms  that  play  important  roles  in  describing  the  nonlinear  evolution  of  disturbances. 
Such  an  interacting  triple  is  possible  whenever  three  waves  with  wavenumbers  k,  and 
frequencies  Qj  {j  =  1,2,3)  such  that  ki  =  k2  +  ka  and  fli  =  fl2  +  fla  can  coexist. 
This  type  of  mechanism  is  the  subject  of  our  current  investigation,  but  an  excellent 
account  of  the  variety  of  interaction  phenomena  which  are  possible  within  the  scenario 
of  fluid  mechanics  has  been  given  by  Craik  (1). 

The  crossflow  vortex  instabihty  structure  considered  here  can  only  occur  in  three- 
dimensional  boundary  layers  and  was  first  examined  both  theoretically  and  experimen¬ 
tally  by  Gregory,  Stuart  &  Walker  (2)  in  the  context  of  rotating  disc  flow.  They  showed 
with  the  china-clay  technique  that  this  instability  takes  the  form  of  a  regularly  spaced 
pattern  of  equiangular  spiral  vortices  which  are  stationary  relative  to  the  disc.  In  this 
paper  GSW  showed,  using  inviscid  theory,  that  the  stationary  vortex  mechanism  is  as¬ 
sociated  with  an  inflexional  profile  in  which  the  inflexion  point  coincides  with  a  point 
of  zero  velocity  somewhere  in  the  flow.  Their  calculation  of  the  number  of  vortices 
expected  to  be  seen  in  the  flow  above  a  rotating  disc  gave  a  number  approximately 
four  times  greater  than  the  observed  value  of  about  30,  but  the  angle  of  13®  between 
the  axes  of  the  vortices  and  the  radius  vector  predicted  was  in  good  agreement  with 
that  seen  in  experiments.  Brown  (3),  and  later  Cebeci  &  Stewartson  (4)  solved  the 
Orr-Sommerfeld  equation  numerically  and  obtained  values  for  the  critical  Reynolds 


number  for  transition  much  lower  than  that  observed  by  Gregory  et.  al.  and  used 
by  Stuart  in  his  calculation  for  the  number  of  vortices.  Malik  (6)  found  the  neu¬ 
tral  stability  curve  for  these  stationary  disturbances  and  he  showed  the  existence  of 
a  second  stationary  mode  (this  one  however  being  viscous)  which  corresponds  to  zero 
wall  shear  stress  of  the  crossflow  velocity  profile.  The  asymptotic  structures  (at  high 
Reynolds  numbers)  of  both  this  viscous  mode  and  of  the  inviscid  mode  of  GSW  were 
obtained  by  Hall  (6).  The  linear  analysis  presented  in  that  paper  has  been  extended 
by  MacKerrell  (7)  and  Gajjar  (8)  to  include  nonlinear  effects. 

The  aforementioned  investigations  have  concentrated  on  stationary  modes.  How¬ 
ever,  Faller  &  Kaylor  (9)  and  Lilly  (10),  in  their  studies  of  periodic  perturbations 
imposed  on  Ekman  boundary  layer  flows,  mentioned  the  existence  of  an  instability 
with  non-zero  phase  speed.  Federov  et.  al.  (11)  showed  that  in  rotating  disc  flow 
a  similar  instability  occurs  much  earlier  than  the  stationary  mode  observed  in  the 
experiments  of  Gregory  et.  al..  This  instability  also  appeared  as  a  pattern  of  spiral 
vortices,  but  which  was  moving  relative  to  the  disc.  The  number  of  these  vortices  was 
observed  to  be  between  14  and  16,  and  had  axes  inclined  at  about  20®  to  the  radius 
vector.  Recently,  Bassom  &  Gajjar  (12)  have  conducted  a  theoretical  analysis  of  some 
aspects  of  non-stationary  vortices  in  rotating  disc  flow  and,  in  particular,  investigated 
the  properties  of  neutrally  stable  nonlinear  modes.  Our  present  concern  is  with  the 
time  dependent  version  of  the  high  Reynolds  number  inviscid  crossflow  vortex  whose 
flow  structure  wjis  elucidated  by  Hall  (6).  The  modifications  required  to  his  anal¬ 
ysis  to  incorporate  the  effects  of  unsteadiness  are  relatively  straightforward  and  are 
considered  in  the  following  section. 

There  has  been  much  work  in  relation  to  the  important  problem  of  the  interaction 
of  Tollmien-Schlichting  disturbances  with  Gortler  vortices  (an  instability  associated 
with  flows  over  curved  surfaces).  See,  for  example,  (13-17)  and  the  references  therein. 
Also,  some  study  has  been  made  of  the  interaction  between  Tollmien-Schlichting  waves 
and  steady  crossflow  vortices,  (18-20),  and  of  the  triad  interaction  between  three  un¬ 
steady  crossflow  vortices  under  consideration  here.  In  particular,  we  refer  to  the  work 
of  Reed  (18,  19)  and  El-Hady  (21),  who  found  that  the  growth  of  each  of  the  insta¬ 
bility  modes  could  be  substantially  influenced  by  their  mutual  interaction  mechanisms. 
Additionally,  non-parallel  effects  were  included  and  it  was  concluded  that  the  evolu¬ 
tion  of  the  vortices  was  critically  dependent  upon  the  initial  amplitudes  of  the  modes. 
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However,  these  accounts  ignored  the  crucial  fact  that  for  a  completely  rational  descrip¬ 
tion  of  the  importance  of  the  interaction  process  in  the  development  of  the  disturbances 
each  of  the  instability  modes  involved  should  be  neutrally  stable  at  leading  order  in 
their  own  rights.  Otherwise,  the  vortex  growth  due  to  the  interaction  mechanism  is 
no  larger  than  would  be  experienced  by  the  modes  in  the  absence  of  the  interaction, 
and  the  true  importance  of  the  mechanism  is  extremely  difficult,  if  not  impossible,  to 
accurately  assess.  For  this  reason,  our  present  aim  is  to  first  illustrate  that  a  suitable 
triad  of  travelling  crossflow  vortices  can  exist  in  a  three-dimensional  boundary  layer; 
crucially,  these  modes  are  chosen  so  as  to  ensure  that  each  is  individually  neutrally 
stable  at  leading  order. 

Once  the  existence  of  the  desired  triad  is  established  we  move  on  to  consider  the 
problem  when  all  three  modes  are  simultaneously  present.  In  the  case  of  a  single, 
inflnitesimally  small  mode  the  asymptotic  flow  structure  consists  of  a  main  layer  (co¬ 
incident  with  the  majority  of  the  boundary  layer)  augmented  by  a  thin,  viscous  wall 
layer  next  to  the  surface  of  the  disc  and  by  at  least  one  linear  critical  layer  present 
somewhere  in  the  flow.  With  all  modes  of  the  triad  present,  the  corresponding  flow 
structure  becomes  much  more  complicated,  for  now  we  have  at  least  three  critical  lay¬ 
ers  (one  fo;  each  mode)  and  sometimes  more.  It  transpires  that  in  order  to  obtain  a 
rational  description  of  the  evolution  of  the  triad  each  of  the  critical  layers  needs  to  be 
carefully  analysed,  for  each  contributes  non-zero  terms  to  each  of  the  final,  coupled 
amplitude  equations  for  the  three  modes. 

An  important  characteristic  of  our  analysis  is  that  each  vortex  is  taken  to  be  suffi¬ 
ciently  small  and  to  develop  on  appropriately  slow  length  and  timescales  such  that  all 
the  critical  layers  remain  linear  in  character.  We  emphasize  that  although  the  equa¬ 
tions  obtained  in  this  paper  are  specifically  oriented  towards  the  rotating  disc  flow,  for 
the  cases  of  other  basic  velocity  profiles  the  derivation  procedure  would  be  essentially 
identical.  The  method  adopted  also  allows  nonparallel  effects  to  be  accounted  for  with¬ 
out  formal  difficulty.  The  linked  triad  of  evolution  equations,  somewhat  unfortunately, 
suffer  from  the  handicap  that  they  are  very  complicated  and  require  much  delicate  non¬ 
trivial  numerical  work  to  establish  the  interaction  coefficients  for  any  particular  flow. 
Here  we  perform  some  of  the  calculations  for  the  rotating  disc  situation  to  demonstrate 
that  the  derivation  of  the  evolution  equations  is  practicable.  The  equations  so  derived 
have  complex-valued  coefficients  and  are  not  dissimilar  in  form  to  those  derived  by 
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Smith  &:  Stewart  (22)  who  were  concerned  with  high  frequency  resonant  triads  within 
boundary  layers.  We  consider  that  the  principal  result  of  the  current  work  is  that  in 
three-dimensional  boundary  layers  resonant  triads  of  non-stationary  crossflow  vortices 
do  exist.  Further,  given  such  a  triad  of  modes  which  are  neutrally  stable  at  leading  or¬ 
der,  we  show  that  we  can  develop  a  fully  rational  asymptotic  account  of  the  interaction 
processes. 

The  procedure  for  the  remainder  of  the  paper  is  as  follows.  In  section  2  we 
briefly  consider  the  asymptotic  structure  of  a  Unear,  travelling  crossflow  vortex  in  a 
three-dimensional  boundary  layer,  making  specific  reference  to  the  case  of  rotating 
disc  flow  as  an  example.  Numerical  work  in  section  3  estabUshes  that  the  desired  form 
of  a  resonaiit  triad  of  modes  exists  for  this  special  basic  flow  and,  armed  with  this 
information,  we  make  a  careful  study  of  the  interaction  mechanism  in  section  4  to 
enable  us  to  derive  a  coupled  triple  of  evolution  equations  for  the  modes.  Finally,  we 
perform  some  numerical  work  to  estabUsh  the  interaction  coefficients,  discuss  a  few 
solutions  of  the  relevant  evolution  equations  and  we  conclude  with  some  discussion 
and  suggestions  for  further  work. 


§2.  The  basic  equations  and  flow  structure  for  a  linear  crossflow  vortex. 

We  have  already  stated  that  we  intend  to  make  an  investigation  of  the  interaction 
processes  which  take  place  between  the  three  members  of  a  resonant  triad  of  crossflow 
vortices  and  a  natural  first  step  is  to  consider  the  disturbance  structure  for  each  indi¬ 
vidual  mode  in  isolation.  Here  we  concentrate  on  the  rotating  disc  flow  problem  for 
definiteness,  although  the  work  can  be  trivially  modified  to  account  for  other  three- 
dimensional  basic  flows;  for  example  the  asymptotic  suction  profile  or  the  attachment 
line  flow,  see  Hall  Si  MacKerrell  (23).  In  the  particular  problem  at  hand,  we  consider 
the  case  in  which  the  disc  rotates  about  the  2-axis  with  angular  velocity  fl.  Relative 
to  cylindrical  polar  co-ordinates  (r,&,  z)  which  rotate  with  the  disc  and  in  which  r  and 
2  have  been  made  dimensionless  with  some  reference  lengthscale  L,  the  continuity  and 
Navier-Stokes  equations  for  an  incompressible  fluid  in  the  region  2  >  0  are 

V.u  =  0,  (2.1a) 

dVL  1  o 

—  (u.V)  u  -I-  2  (fi  A  u)  +  n  A  (n  A  r)  =  — Vp  -|-  t^V^u.  (2.16) 
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Here  u  is  the  velocity  vector,  r  is  the  co-ordinate  vector,  p  is  the  fluid  density,  p  is  the 
fluid  pressure  and  u  the  kinematic  viscosity.  The  Reynolds  number  R  for  the  flow  is 
given  by  i?  =  ^l?‘lv  and  is  taken  to  be  large  throughout  the  following. 

It  is  found  convenient  to  define  the  small  parameter  e  =  R~  • .  As  the  axes  rotate 
with  the  disc,  the  basic  flow  is  given  by  the  Von-KArmdn  solution 

u  =  Ub  =  Tn  (r’0(0,r{5(^),c^u)(^))  ,  p  =  pL‘^n^€^p{i),  (2.2) 

where  z  =  and  U,v,  w  satisfy  the  equations 

-  {1  +  d)'^  +  n  w  =  u\  2q(1  +  P)  +  o'ffi  =  t)",  (2.3a, 6) 

w'+2u  =  Q.  (2.3c) 

Here  primes  denote  differentiations  with  respect  to  ^  and  the  appropriate  boundary 
conditions  are 


U  =  P  =  tjj  —  0  on  4  =  0;  U  — »  0,  V  — *  —1  as  4  — ^  oo- 


The  structure  of  an  infinitesimal  unsteady  crossflow  vortex  imposed  on  the  above 
basic  flow  follows  closely  that  given  by  Hall  (6)  for  the  steady  case.  The  mode  has 
wavelengths  scaled  on  the  boundary  layer  thickness  and  an  0(1)  wavespeed  and  so  we 
consider  perturbation  quantities  with  r  and  9  dependence  given  by  E,  defined  by 


E  =  exp 


(/' 


a{r,€)dr  +  90{e) 


(2.4a) 


We  shall  restrict  our  attention  to  examining  disturbances  in  a  neighbourhood  of 
some  point  r  =  and  expand  the  wavenumbers  and  frequency  in  the  forms 


a  =  ao  +  eai  +  . . . ,  /?  =  /?o  +  f/?i  +  . . . ,  =  flo  +  efli  + -  (2.46) 

In  general  these  quantities  will  be  complex  but  here  we  shall  be  concentrating  on 
neutrally  stable  disturbances  (at  least  to  leading  order)  and  so  we  take  ao,  /^o  and  flo 
to  be  real.  The  disturbance  structure  in  the  direction  is  sketched  in  figure  (1).  There 
is  an  inviscid  layer  (/)  of  thickness  O(e^)  (the  same  depth  as  the  boundary  layer),  and 
to  satisfy  no-slip  conditions  at  the  wall  a  viscous  layer  {II)  of  thickness  O  (e^^^)  must 
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be  present.  If  the  basic  flow  is  perturbed  so  that  the  total  non-dimensional  flow  field 
is  given  by 

(u,  V,  w,p)  —  {rU,r'D,  +  6  [(U,V,  W,P)  E  c.c.]  +  . . . ,  (2.5a) 

where  6  1,  then  in  zone  (/)  the  disturbance  quantities  U,V,W,P  (which  are  func¬ 

tions  of  ^  alone)  expand  in  the  forms 

{U,V,W,P)  =  iUQ  +  €Ui+  ...,Vq  +  €Vi  +  ...,Wo  +  €Wi  +  ...,Po  +  €Pi  +  ...). 

(2.56) 

Substituting  expressions  (2.5)  in  the  continuity  and  the  linearised  Navier-Stokes 
equations  yields  at  leading  orders 

i  ( aoUo  -I-  ^ Vb)  +Wi  =  0,  (2.6a) 

V  r„  / 

i  {aornU  +  /3o0  -  IIo)  Uo+  r^U  Wq  =  -iaoPo,  (2.66) 

I  i3c\ 

i  (aor„G  +  /?oG  —  fIo)^0  +  ^0  = - -POt  (2.6c) 

fn 

and 

*  («0fnQ  +  ~  flo)  ^0  =  — Pq-  (2.6d) 

If  we  define 

Uj  =  a,r„s  ■¥  0j’D-  n,-, 
then  we  obtain  the  classical  Rayleigh  equation 

{Wq  -  =  V'q  Wo,  (2.7a) 

where  79  =  oq  +  The  required  boundary  conditions  on  (2.7a)  are,  firstly, 

|W^0(oo)|  <  00,  I^o(0+)  =  0,  (2.76,  c) 

to  avoid  exponential  growth  at  the  edge  of  the  boundary  layer  and  to  give  the  typical 
inviscid  condition  of  tangential  flow  as  the  wall  is  approached.  Technically,  to  reduce 
this  slip  velocity  to  zero  at  the  surface  of  the  disc,  analysis  of  the  viscous  layer  {II)  is 
required  but  such  analysis  does  not  affect  the  leading  order  quantities  oq,  00  and  flo 
which  are  of  the  primary  interest  here. 
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The  Rayleigh  equation  (2.7a)  additionally  has  singularities  at  any  positions  where 
Uq  =  0,  and,  if  this  occurs,  the  presence  of  an  internal  critical  layer  (at  that  position) 
is  suggested.  The  details  of  the  asymptotic  analysis  of  such  a  (linear)  critical  layer 
follow  well-known  lines  (see  for  example  Drazin  k.  Reid  (24))  and  so  we  merely  sketch 
a  few  of  the  relevant  results  here.  If  in  the  vicinity  of  the  critical  layer,  at  ^  ^  say, 

we  have  the  Taylor  expansion 

Oi=i.Si^U-0'  ,  (2.8a) 

-1 
f  =1 

where  we  have  defined  |j^-,  then  the  solution  of  (2.7a)  as 

i  is  given  by 

00  oo 

O'  +  T>  Grii- O' lriU-0  ,  (2.8b) 

T-0  T=1 

where  are  constants  upon  which  the  remaining  coefficients  (n  =  2,3,...), 

G„  (n  =  1,2,...)  depend.  In  particular,  we  have  Gi  =  goBo2/^oi-  The  solution 
(2.86)  implies  that  Bq  and  Vq  in  (2.5)  have  logarithmic  singularities  as  ^  — *■  (+  and 
to  smooth  these  we  need  to  invoke  viscous  effects  in  a  thin  layer  (III  in  figure  1) 
surrounding  ^  —  0  Balancing  inertial  and  diffusive  terms  suggests  that  this  zone  is 
of  thickness  G(e^)  and  if  we  define  ^  ^  +  eip  here,  where  ip  =  0(1),  the  disturbance 

quantities  develop  as 

U  =  €~  ^Uq  +  In  eOi  +  +  eln  +  •  •  • ,  (2.9a) 


V  =  e"  +  In  eVi  +  V2  +  eln  +  eV^  +  . . . ,  (2.96) 

W  =  Wq  +  e In  eWi  +  in  (W^  +  +  . . . ,  (2.9c) 

P  =  Po  +  elnePi  +  eP2  +  c^lnePs  +  e^P^  +  —  (2.9d) 

Briefly,  we  find  on  substituting  (2.9)  into  the  linearised  Navier  Stokes  equations 


that 

Uo  = 


\Boi\ 


ao - - - I  Pq  /  exp 

Boi  /  •'0 


irasn(Soi)  (’/>  +  |^)]  (-JJ^)  dr. 


3|Boil' 


(2.10a) 
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where  Pq  =  Po(0-  This  solution,  as  |0|  — *  oo,  matches  with  the  solution  for  Uq  in 
(2.6)  as  ^  Across  the  critical  layer  we  also  find  that 


Wo=- 

hoPo 

001  ’ 

(2.106) 

ji,  50002 

Bo, 

(-1=) 

,  (2.10c) 

and 

d^W2 

dip^ 

_  i  .00250 
|0oi| 

f  exp  —iTsgn{BQi) 
0 

('•if) 

If  we  integrate  (2.10d)  once  and  then  apply  the  asymptotic  results  given  by  Haber- 
man  (25)  concerning  the  behaviour  of  integrals  of  the  general  type 


for  large  argument  ip,  we  may  show  that  the  derivative  dW2/dip  suffers  a  discontinuity 
across  the  critical  layer  and  in  particular  that 


r  "  1 

dW2 

.  dip  . 


OO 


-  00 


inBo2go 
l^oil  ■ 


(2.11) 


The  implication  of  the  solutions  (2.10),  (2.11)  is  that  the  Frobenius  solution  (2.86) 
for  the  eigenfunction  Wq  as  ^  — »  (+  continues  to  hold  in  the  limit  ^  so 

long  as  the  replacement  of  ln(^  —  ^)  by  (ln(|^  —  ^)  —  iirsgn{Boi))  in  (2.86)  is  made. 
Equivalently,  we  have  that  as  ^  — >  (  —  , 


OO  OO 

H  +  IZ  G'r  [ln(f  -  0  -  iTrsgn{Boi)]  ,  (2.12) 

r  =0  T=1 


just  below  the  critical  layer.  This  is  the  well  known  phase  shift  of  ‘eV’  across  a  linear 
critical  layer. 

We  now  have  a  complete  specification  of  the  problem  to  determine  the  leading 
order  neutral  wavenumbers  qq,  /?o  at  some  given  frequency  flo.  The  Rayleigh  equation 
(2.7o)  needs  to  be  solved  subject  to  the  boundary  conditions  (2.76,  c)  and  where  the 
replacement  given  by  (2.12)  in  the  local  series  solution  (2.86)  across  the  critical  layers 
is  applied.  The  analysis  presented  here  is  specific  to  the  case  of  rotating  disc  flow, 
but  may  easily  be  modified  to  account  for  a  variety  of  other  three-dimensional  basic 
velocity  profiles. 
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§3.  The  isolation  of  resonant  triads. 

At  the  outset  our  aim  was  to  verify  that  resonant  triads  (with  the  property  that 
each  member  of  the  triad  is  neutrally  stable  at  leading  order)  are  possible  within 
particular  three-dimensional  boundary  layer  flows.  In  the  context  of  rotating  disc 
flow,  we  therefore  wish  to  seek  ordered  triples  {(ooj » /^Oj » )  \j  =  15  2,3}  such  that 

^^01  =  Q!02  +  Q!03.  PO1  =  0O2  +  0O3^  f^Ol  =  1^02  +  flo3>  (3.1) 

and  with  the  characteristic  that  each  triple  is  an  eigenvalue  set  of  the  system  (2.7), 
(2.8)  &  (2.12).  To  this  end  a  numerical  approach  was  adopted. 

A  flrst  concern  is  to  And  the  location(s)  of  the  critical  layer(s).  These  occur 
wherever  Oq  =  0;  i.e.  wherever 


aor„  u  + (3.2a) 

or,  if  we  define  the  leading  order  wavespeed  cq  =  wherever 

17,  =  G  + /xC  =  CO,  (3.26) 

with  (j,  =  /5o/(r„ao).  Hence,  the  ‘effective’  basic  velocity  profile  U,  is  a  function  of  n 
and  this  flow  is  sketched  in  figure  (2).  It  is  known  that  q'(0)  =  0.510,  G  (0)  =  —0.616 
and  then  as  the  parameter  /x  varies  the  effective  basic  velocity  acquires  a  number  of 
distinct  possible  forms.  The  boundary  conditions  in  (2.3)  suggest  that  U,  — >  — /x  as 
^  — >  oo  and  for  /x  >  /x^  (w  0.8284)  the  profile  is  a  monotonic  decreasing  function. 
For  0  <  n  <  fic,  U,  vanishes  somewhere  other  than  at  the  origin,  and  indeed  when 
H  «  0.235,  the  point  of  inflexion  of  the  basic  flow  U,  coincides  with  the  point  at 
which  it  vanishes.  In  this  case  the  boundary  layer  becomes  susceptible  to  the  inviscid 
stationary  vortex  found  by  GSW.  For  /x  <  0,  17,  is  always  a  positive  function  and  for 
/X  <«  —0.6  it  is  monotonic  as  well.  We  see  from  figure  (2)  that,  depending  upon  the 
values  of  /x  and  cq  in  (3.26),  v’e  may  have  either  one  or  two  critical  layers. 

The  numerical  scheme  used  in  practice  requires  that  a  value  of  ^  be  preselected 
and  then  ao,  /3o  and  Qq  are  determined  such  that  /x  =  /?o/(r„ao).  It  is  easily  shown 
from  (2.7a)  that  as  ^  — >  oo  the  eigenfunction  Wq  is  proportional  to  exp(— 70^). 
A  suitably  large  value  of  say  ^00 ,  was  chosen  and  the  asymptotic  behaviour  of 
the  eigensolution  was  assumed  to  hold  there.  With  these  starting  conditions,  the 
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Rayleigh  equation  (2.7a)  was  integrated  with  a  fourth  order  Runge-Kutta  scheme  to 
a  point  a  small  distance,  say  A,  above  the  first  critical  layer  encountered.  There  the 
series  solution  (2.86)  was  assumed  to  hold  and  qq  and  g\  were  chosen  so  that  the 
necessary  matching  conditions  held.  Then,  using  the  replacement  suggested  by  (2.12), 
the  solution  just  below  the  critical  layer  was  established  which  could  then  be  used  to 
integrate  (2.7a)  until  either  the  next  critical  layer  position  was  reached  (in  which  case 
the  process  described  above  was  repeated)  or  until  ^  =  0.  In  general,  the  predicted 
value  of  iyo(^  =  0)  is  complex  and  non-zero  and  so  a  Newton  iteration  on  the  real 
values  oq  cq  was  employed  to  ensure  that  the  required  condition  (2.7c)  was  indeed 
met.  Given  the  eigenvalues  ao  and  cq  for  some  chosen  value  /x  it  is  then  a  simple 
task  to  retrieve  the  leading  order  azimuthal  wavenumber  /?o  and  the  frequency  flo- 
This  numerical  work  was  checked  for  independence  from  the  chosen  value  ^oo  ,  the  step 
length  used  in  the  Runge-Kutta  integration  and  of  the  distance  A  from  the  critical 
layer(s)  at  which  the  series  solutions  were  applied. 

The  solutions  of  the  eigenproblem  for  various  values  of  the  parameter  /x  (which 
is  tan~^<^,  where  (p  is  the  angle  which  the  considered  disturbances  make  with  the 
outward  radial  direction  on  the  disc)  are  shown  in  figure  (3),  where  we  illustrate  the 
dependence  of  the  leading  order  radial  wavenumber  ao  and  the  frequency  flo  on  the 
wavenumber  /3o-  There  are  several  noteworthy  features  of  the  solutions.  Firstly,  there 
are  two  distinct  branches  in  solution  space  which  are  independent  of  each  other  and  do 
not  coalesce.  On  one  branch  (labelled  /  on  figure  (3))  we  see  that  there  is  a  non-trivial 
solution  of  the  eigensystem  with  zero  frequency  and  this  corresponds  to  the  inviscid 
stationary  mode  of  (2).  Here  qq  %  1.1,  /?o  ^  0.26r„  and  /x  «  0.235,  and  then  the  basic 
profile  U,  has  a  point  of  inflexion  coincident  with  the  point  at  which  it  vanishes.  As 
the  parameter  /x  increases  we  find  that  the  frequency  becomes  negative  and  using  the 
sketches  of  the  basic  flow  in  figure  (2)  we  see  that  we  are  reduced  to  an  eigenproblem 
with  a  single  critical  point.  Appealing  to  the  ideas  of  Stuart  in  (2),  in  the  case  of 
only  one  critical  layer  the  neutral  wavespeed  cq  is  determined  by  the  condition  that 
the  function  U,  in  (3.26)  must  have  a  point  of  inflexion  where  C/,  =  cq. 

For  a  waveangle  parameter  fj,  less  than  the  inviscid  stationary  mode  value  of  0.235 
we  now  have  two  critical  layers  in  the  flow  and  the  neutral  frequency  is  positive.  As 
fi  — *  0  the  wavenumbers  and  frequency  become  small.  Finally  for  this  branch,  as 
/X  — >  0  or  /X  — >  oo  then  the  wavenumbers  in  the  radial  and  azimuthal  directions 
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become  of  different  orders  of  magnitude  and  the  symmetry  of  the  governing  equations 
(2.6)  is  lost.  The  structures  which  would  govern  these  limiting  cases  are  then  not 
covered  by  the  work  of  section  2  and  it  would  be  interesting  to  identify  and  analyse 
these  structures. 


On  the  second  branch  of  solutions  (labelled  II  on  figure  3)  it  is  immediately  seen 
that  only  positive  frequencies  are  found.  In  particular,  we  remark  that  as  /i  increases 
and  approaches  Hc  ^  0.187  (corresponding  to  a  disturbance  waveangle  of  10.6°)  the 
wavenumbers  and  frequency  become  small.  This  suggests  the  existence  of  a  longer 
wavelength  disturbance  with  a  small  wavespeed  and  this  structure  corresponds  to  that 
of  the  ‘upper-branch’  mode  investigated  in  (12).  In  this  paper,  small  perturbations 
with  wavenumbers  O(iJ^)  and  wavespeeds  0{R~'^)  were  shown  to  be  possible  at 
waveangles  between  approximately  10.6°  and  39.6°.  As  /r  decreases  from  fic  in  our 
current  problem  and  becomes  negative  the  convergence  of  the  numerical  scheme  em¬ 
ployed  becomes  increasingly  difficult  to  achieve  and  only  the  neutral  solutions  which 
can  be  reliably  located  by  our  method  have  been  shown  in  figure  (3).  Indeed,  for 
sufficiently  negative  /x  the  basic  profile  is  monotonically  increasing  and  has  no  point 
of  inflexion.  At  these  values  of  parameter,  no  neutral  solutions  of  the  sought  type  can 
exist  and  so  it  suggests  the  existence  of  some  <  0  at  which  our  solutions  cut  out.  It 
is  seen  from  figure  (2)  that  as  fi  decreases  from  zero  the  neutral  disturbances  have  two 
critical  layers  up  to  the  point  where  the  basic  velocity  is  no  longer  non-monotonic. 
Then  the  two  critical  layers  must  merge,  cease  to  remain  asymptotically  distinct  and 
a  fresh  analysis  is  required.  Since  the  principal  thrust  of  the  present  work  is  to  study 
resonant  triads  we  do  not  explore  this  aspect  any  further  here.  However,  we  remark 
that  these  considerations  show  that  neutrally  stable  linearised  disturbances  are  pos¬ 
sible  at  propagation  angles  >  tan"^(/x.).  There  is  therefore,  a  lower  bound  for  this 
angle  but  no  formal  upper  bound  since  from  figure  (3)  we  can  see  that  neutral  modes 
exist  for  /x,  <  /x  <  oo.  This  contrasts  with  the  results  of  (12)  where  it  was  found  that 
for  linear  disturbances  the  possible  waveangles  associated  with  upper  branch  neutral 
modes  were  bounded  both  above  and  below. 


Given  the  eigenvalues  of  the  problem  specified  by  (2.7),  (2.8)  &  (2.12)  we  sought 
resonant  triad  triples  (ooj  ^  (^Oj  >  waves  which  are  each  neutrally  stable  at  leading 


11 


order  and  which  satisfy  requirements  (3.1).  We  found  two  such  sets,  in  particular 

aoi  =0.8514,  ^01  =  0.2648r„,  Jloi  =  -3.137  x  IQ-^, 

Q02  =  0.7343,  ;3o2  =  0.0731r„,  ^02  =  1-12  x  10"  2,  (3.3a) 

Q03  =0.1171,  /3o3  =  0.1917r„,  flo3  =  -4.25  x  lO-^, 

and 

aoi  =0.5847,  /?oi  =  0.2764r„,  ^oi  = -6.29  x  lO'^, 

002  =  0.36  05,  /3o2  =  0.2699r„,  flo2  =  -7.16  x  10" 2,  (3.36) 

003  =0.2242,  /3o3  =  6.46  X  10~^r„,  ^03  =  8.76  x  10'^. 

We  do  not  claim  that  these  are  the  only  possible  triads  but  the  important  result 
here  is  that  resonant  triples  with  the  sought  neutral  stability  properties  do  indeed  exist 
within  a  three-dimensional  boundary  layer.  The  type  of  analytical  and  numerical  work 
presented  here  can  be  applied  to  other  basic  velocity  profiles.  Indeed  very  similar 
results  as  to  the  forms  of  the  neutral  stability  curves  in  oo//?o  and  flo/j^o  spaces 
appear  for  the  attachment  line  profile  and,  of  course,  given  the  existence  of  triads 
for  a  simple  rotating  disc  flow  we  have  no  reason  to  doubt  that  a  similar  conclusion 
holds  for  other  much  more  complicated  and  more  physically  realistic  three-dimensional 
boundary  layers. 

This  verification  of  the  existence  of  resonant  triads  leads  us  on  to  an  integral 
part  of  the  current  work-  that  of  considering  the  disturbance  structure  of  the  three 
members  of  the  triad  when  they  are  permitted  to  interact  and  volve  on  a  suitably 
slow  timescale. 


§4.  An  asymptotic  description  of  the  interaction  mechanisms  in  a  resonant 
triad  of  crossflow  vortices. 

We  have  illustrated  that  within  the  context  of  three-dimensional  boundary  layer 
flows  resonant  triads  comprising  of  crossflow  vortices  can  arise.  We  now  investigate 
an  asymptotic  description  of  how  the  members  of  the  triad  may  evolve  in  time  and 
space  once  the  interaction  processes  are  accounted  for.  To  ensure  that  the  structures 
of  the  critical  layers  in  the  flow  remain  as  simple  as  possible  for  this  study  we  assume 
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that  the  amplitude  of  each  vortex  is  0(S)  where  (5  <  1  is  suiSciently  small  such  that 
all  the  the  critical  layers  are  linear  in  character.  Then,  through  the  nonlinear  terms 
in  the  Navier  Stokes  equations,  it  is  found  that  each  pair  of  modes  combine  to  drive 
the  third  mode  via  terms  of  size  O(S^)  and  this  in  turn  suggests  that  the  interaction 
mechanisms  force  the  vortices  to  evolve  on  slow  length  and  time  scales  of  size  0(6). 

If  again  we  relate  our  working  to  the  rotating  disc  problem  as  a  concrete  example, 
suppose  that  we  wish  to  describe  the  interaction  of  three  hnear  vortices  with  leading 
order  wavenumbers  ao, ,  /3oj  and  of  frequencies  Qoj  (j  =  1(2,3)  such  that  the  resonance 
conditions  (3.1)  hold.  Then  we  define 


aj  (r,  €)dr  +  0/8j  (e)  -  tflj 


(4.1a) 


where 


«;■  =  O!0j  +  eaiy  +  . . . ,  j3j  =  /3oj  +  efSij  +  ...,  Q,  =  flo;  +  +  . . . ,  (4.16) 


and  El  =  E2EZ.  We  allow  the  disturbance  vortices  to  develop  on  slow  lengthscales 
ri  =  6  (r  —  r„),  di  =  66  and  timescale  ti  =  6t,  so  that  in  the  governing  Navier  Stokes 
equations  we  formally  use  the  methods  of  multiple  scales  to  replace  the  r,  9  and  t 
derivatives  according  to  the  transformations 


^ ^  _a  _a_ _ ^  d 

dr  dr  dri  dO  d9  d9i'  ,  . 

a_ _ ^d_ 

dt  dt  ^  dti 

As  previously,  we  choose  to  examine  vortices  in  the  neighbourhood  of  some  point 
r  =  r„  and  then,  in  the  main  inviscid  region  (I)  where  z  =  we  are  led  to  consider 
perturbation  quantities  of  the  forms 


U  =  6 


(^00;  +  •  •  •)  -A-jEj  +  C.C. 

L;-l 


Lj-1 


W  =  6 


1.  J 


1 


(Wooy  +  • . .)  ^}Ej  -|-  c.c.j  +  <5^ 


L;=l 


YL  ^10;  •  +  .  .  .  j  +  .  .  .  , 

(•i.3a) 

;=1  J 

3  1 

Y1  ^10;  +  .  .  .  1  +  .  .  .  , 

(4.34) 

,;=1  J 

1  ^  1 

E  Wio^E^+...\  +..., 

[j=l  J 

(4.3c) 
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In  (4.3)  the  scaled  amplitudes  Aj  are  functions  of  the  slow  scales  ri,  0i  and  ti 
and  Uooj ,  VbO;  >  •  •  •  >  >  •  •  •  are  functions  of  ^  alone.  Substituting  (4.2)  and  (4.3)  into 

the  Navier  Stokes  equations  and  comparing  coefficients  of  Ej  we  find,  as  in  section  2, 
that  at  leading  orders  FToOj  satisfies  the  Rayleigh  equation 


Uok  {WoQ,  -  7?*  W^OO*)  -  Uq,  W^oo*  =0,  A:  =  1, 2,  3,  (4.4) 


where  =  aji,rnU+Pjkd—Q,jk  and  7?^  =  +  (/?/t /»'^)-  The  functions  Uooj  and  Vboj 

may  be  found  in  terms  of  VFqOj  by  suitable  manipulations  of  the  governing  equations, 
which  are  precisely  (2.6)  with  ao,/?0)f^O  replaced  by  ao; , /^O;  >  ^^0;  and  Uo,Vo,Wq,Po 
by  Uooj ,  Vooj ,  kFoOj ,  • 

^From  here  onward  we  concentrate  on  deriving  the  evolution  equation  for  the 
vortex  proportional  to  Ei.  Given  the  obvious  symmetry  in  the  problem,  together 
with  the  generalised  way  in  which  the  amplitude  equations  will  be  derived,  once  the 
governing  equation  has  been  obtained  for  one  mode  it  is  formally  straightforward 
to  obtain  the  corresponding  equations  for  the  other  two  modes.  At  0(6^)  in  the 
Navier  Stokes  equations  we  find  that  the  function  Wioi  in  (4.3c)  satisfies  the  non- 
homogeneous  Rayleigh  equation 

C^Oi  (^101  -  701^101)  -  = 

-  iWooi  “  2iUoiWooi  (aoi^  +  A  A2A3, 

(4.5a) 

where 


I*  =  i  (/9o30!01  +  ao2l3oi)  +  (^02001  +  ao3/?Ol)  +  aQiUoo2l 

rn 

^^Vbo2^oo3  +  W^oo3  (0:01^^002  +  ^^Vbo2)  +  1^002  (Q;oii^oo3  +  ^^^003) 

J  L  ^  r„  /  ^  r„  / 

—  *701  *  (  002^^003^^002  +  003^003^^002  +  Vo03^002  +  ^^Vbo2W^003) 

>■  \  r„  / 

+  W003^002  +  ^002^003]  ’ 


+  Ooii7o02t^003+ 


(4.56) 
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and 


^  f  9  d  d  \  ,  ^ 

The  appropriate  boundary  conditions  on  the  function  Wqi  are  very  similar  to 
those  for  VFoOl  in  (4-4).  In  particular,  to  avoid  unbounded  growth  as  ^  — >  oo  requires 
that  |iyioi(cx))l  <  oo  and  for  no  disturbance  at  the  surface  of  the  disc,  1^101(0+)  = 
0.  Additionally,  it  is  clear  that  appropriate  jump  conditions  need  to  be  imposed  at 
the  position  of  the  critical  layer(s).  For  ease  of  notation,  suppose  that  the  mode 
proportional  to  Ej  has  associated  critical  layers  located  at  fy  and  note  that  there 
could  well  be  more  than  one  critical  point  for  each  mode.  Inspection  of  (4.5)  suggests 
that  singularities  in  the  solution  VFioi(0  may  be  expected  at  all  the  critical  layers 
for  j  =  1,2,3,  not  just  at  the  positions 

The  form  of  the  solution  of  (4.5a)  in  the  vicinity  of  the  critical  layers  at  ^y  may 
be  one  of  two  distinct  types.  Close  to  the  position  we  can  easily  show  that  the 
right-hand  side  of  (4.5a)  is  asymptotically  proportional  to  l/(^  —  ^l)  and  from  this  we 
can  deduce  that  for  0  <  ^  —  1, 

IFiOi  =  -Yolog(^  -  ^i)  +  0(1)  +  . . . ,  (4.6a) 


where  Xq  is  a  constant  which  may  be  found  analytically.  Close  to  the  other  critical 
layers,  ^y  {j  =  2,3),  the  behaviour  of  Wioi  is  more  complicated.  Then  we  obtain,  from 
(2.6)  &  (2.8),  that 


T  ,  X\j  X2j 

*  (^-C)"  (e-e,)  ’ 


(4.6i») 


and  then 


W^ioi  ~  ln(^  -  ^y)  +  02y  +  03y(e  -  e";  )ln(e  -  e, )  +  •  •  • ,  (4.6c) 


where  the  constants  Xkj,Dkj  are  given  in  the  Appendix.  The  solutions  (4.6)  are  valid 
for  ^  >  fy  and  clearly  need  to  be  modified  for  the  region  ^  <  ^y .  This  requires  careful 
scrutiny  of  the  critical  layer  properties  in  the  vicinity  of  ^  =  fy  in  much  the  same  way  as 
was  outlined  in  §2.  The  main  difference  now  is  that  the  algebraic  manipulation  needed 
within  the  critical  layer  region  becomes  much  more  complex  due  to  the  presence  of 
all  three  modes.  Rather  than  give  details  of  the  inevitably  lengthy  (but  conceptually 
fairly  straightforward)  critical  layer  analysis,  we  elect  to  merely  state  that  the  result 
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of  this  work  is  as  is  to  be  expected  from  the  classical  theory:  that  is  that  we  obtain  a 
‘itt’  phase  shift  in  the  logarithmic  terms  as  characterised  by  (2.12).  This  then  enables 
us  to  write  down  the  behaviour  of  H^ioi  for  0  <  ^  <C  1  and  thence  to  continue  the 

solution  Wioi  into  the  region  (  <  (j. 

To  derive  the  evolution  equation  for  the  amplitude  Ai  we  follow  the  usual  process 
of  multiplying  (4.5a)  by  14^001  and  integrating  over  the  range  [0,  oo) ,  paying  particular 
attention  to  the  behaviour  of  the  integrals  close  to  the  critical  layers.  Recalling  the 
solutions  (4.6)  and  the  phase  jumps  across  the  critical  layers,  combining  the  integration 
with  the  boundary  conditions  imposed  at  ^  =  0  and  at  ^  leads  to  the  desired 

amplitude  equation.  After  simplification,  and  substituting  the  definition  of  C  given  by 
(4.5c),  we  obtain 


7rr„t2 


7rr„  U  ( 
iBni  1 1  ' 


rj2 

■0012 


r\ 


_  _  -SpiaA 

1-Boiil  -Boil 

^nO"rT/2  tilin' 

-p^W'ooi  - 


finite  part  |5oii| 

{25^12^^001 


1^2  ,  7rRoi2r„QVo2oi 

^^001  +  — 


l-Boill-Boil 


- 


d(  +  2i 


■®011 

»  A 


ooirn  1^001^^^ 


dA\ 

dr\ 


+ 


f  -^012  _  -^013^ 

1-Boiil  ^ 5oii  -Boil' 


TTP 


finite  part  lBoii| 

*^^®-Bo12^001 


W'«ll  + 


001 


oil 


■nBoi2V  H^ooi 

l-Boiil-Boii 


ff 

1 

(  Bqu 

B013]  .„2 

l-Boiil 

-Boil 

•'o  Ool 


dAi 

dti 


(4.7) 


+ 


^  |B0;l| 


A2A3  =  0. 
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In  this  amplitude  equation  we  have  defined 


Bdjk  =  (r„ao,;S^*^(0  + /?0; » 


(4.8) 


and  have  supposed  that 

00 

W^OOl  =  IZ  Qjk  -fi)‘  , 

*=0 

close  to  the  critical  layers  {j  =  2, 3)  at  which  this  function  does  not  have  singularities. 
Otherwise,  the  notation  is  as  previously  defined. 

In  (4.7)  the  expression  ^  denotes  the  Hadamard  (finite)  part  of  the  corresponding 
integral.  Some  of  the  integrals  which  are  contained  in  the  coefficients  of  dAi/dri, 
dAi/ddi  and  dA\/dti  are  singular  at  the  integral  in  the  coefficient  of  ^2^3  is 
singular  at  the  other  critical  points  as  well. 

We  notice  that  if  the  various  integrals  and  summations  were  to  be  evaluated  nu¬ 
merically  and  the  formulae  (4.7)  rearranged  appropriately,  we  could  obtain  an  evolution 
equation  for  the  Ei  mode  which  assumes  the  form; 


dAi 

dt\ 


+  All 


dAi 
dr  I 


+  Ai2 


OAi 


+  Ai3A2i43 


=  0, 


(4.9a) 


where  the  coefficients  Aii,Ai2,  A13  are  some  constants.  By  repeating  the  analysis  out¬ 
lined  above  we  can,  in  an  identical  way,  derive  evolution  equations  for  the  other  two 
modes  of  the  types 


dA2  .  dA2  ,  dA2  ,  . 
dAz  dAz  dAs 


(4.96) 


(4.9c) 


where  an  asterisk  on  a  quantity  denotes  the  complex  conjugate  of  that  quantity. 

Clearly,  the  numerical  work  required  to  actually  evaluate  the  coefficients  in  (4.9) 
for  any  particular  basic  flow  is  nontrivial.  However,  to  illustrate  the  type  of  evolution 
system  which  may  be  obtained  we  considered  the  computational  procedure  relating  to 
the  resonant  triple  given  by  (3.3a).  For  each  member  of  the  triad  we  computed  the 
eigenfunction  Wooj  by  solving  the  Rayleigh  equation  (2.7a)  subject  to  boundary  con¬ 
ditions  (2.76,  c)  and  jump  conditions  (2.12).  The  corresponding  radial  and  azimuthal 
velocity  fields  I/qO;  ,  Vboj  may  then  be  determined  in  terms  of  the  functions  Uq,  and 
Wooj  by  eliminating  the  pressure  from  (2.66,  c).  Since  the  problem  (2.7)  is  hnear  we 
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needed  to  fix  some  normalising  constraint  and  this  was  chosen  such  that  as  ^  — >  oo, 
was  a  real,  positive  function  and  that  max{0  <  ^  <  oo  :  |1^00j  |  =  !}■ 

At  this  stage  the  coefficients  Ayi  and  Xj2  in  (4.9)  were  evaluated.  The  integrals 
involved  in  finding  these  unknowns  are  singular  only  at  and  the  analytic  nature 
of  these  singularities  may  be  deduced  in  a  routine  manner.  To  find  the  finite  parts 
of  the  integrals  across  the  positions  fy ,  we  subtracted  the  singular  behaviour  of  the 
integrand  in  the  neighbourhood  of  and  integrated  the  remaining  (necessarily  well 
behaved)  function  across  .  By  utilising  the  formula  given  by  Lighthill  in  (26)  it 
was  possible  to  deduce  the  Hadamard  part  of  the  subtracted  singularity  function  and 
thence  determine  the  value  of  the  finite  part  of  the  original  integral.  The  method 
worked  well  in  practice  and  the  corresponding  results  are  given  below. 

We  turn  now  to  address  the  much  more  complicated  task  of  determining  the 
interaction  coefficients  Ays  and  again,  as  a  concrete  example,  concentrate  on  the  de¬ 
termination  of  A 13  for  the  method  for  calculating  the  other  interaction  coefficients 
follows  in  a  very  similar  way.  Then  the  function  7,  is  as  defined  in  (4.56)  and  was 
computed  by  using  the  information  previously  derived  on  17oOj  >  k^OO;  and  Wooy  and 
detailed  above.  It  can  be  shown  by  straightforward  analysis  and  using  (4.66)  that  the 
integrand  {LWoqi/Uqi)  has  singularities  at  each  critical  layer  and  to  calculate  this 
integral  a  technique  was  employed  which  is  very  similar  to  that  already  described.  In 
the  vicinity  of  each  critical  layer  the  singular  part  of  the  integrand  was  removed  as 
before,  except  now  we  found  that  the  remaining  function  was  very  erratic  near  the 
critical  point.  This  is  due  to  the  fact  that  the  integrand  can  only  be  determined  nu¬ 
merically  at  the  outset  and  the  process  involves  the  estimation  of  several  derivatives  of 
Uooj,  Vboy  which  is  a  probable  source  of  numerical  error.  On  making  the  usual  checks 
for  numerical  stability,  the  determination  of  the  finite  part  of  the  integral  was  found 
to  be  strongly  influenced  by  the  choice  of  the  intervals  around  the  critical  points  in 
which  the  singular  behaviour  was  removed.  This  is  obviously  unreasonable  and  a  more 
satisfactory  procedure  was  found  in  which  the  integrand  was  written  as  the  sum  of 
the  singular  part  and  an  interpolated  polynomial  in  the  vicinity  of  the  critical  points. 
Even  so,  in  this  case  there  was  found  to  be  some  variation  in  the  determined  value 
of  the  coefficient  A 13,  but  it  has  to  be  remembered  that  this  coefficient  depends  on  a 
calculation  in  which  four  separate  critical  layers  are  traversed  with  a  strong  singularity 
needing  to  be  accounted  for  at  each  critical  layer.  Consequently,  the  difficulties  en- 
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countered  in  obtaining  accurate  values  for  the  interaction  coefficients  in  the  evolution 
equations  (4.9)  are  not  surprising. 

In  this  paper  we  have  managed  to  obtain  approximate  values  for  the  interaction 

^  coefficients.  What  we  have  shown  is  that  determination  of  the  interaction  coefficients 

'  for  the  flow  structure  considered  here  is  possible  in  practice,  and,  with  the  use  of  more 

I  sophisticated  numerical  methods,  there  is  no  reason  to  doubt  that  the  accuracy  of  the 

coefficients  could  be  markedly  improved. 

We  now  describe  some  properties  of  the  solutions  of  the  evolution  equations  and 
conclude  with  some  discussion. 
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§5.  Results  and  Discussion. 


The  numerical  work  described  above  leads  to  evolution  equations  which  assume 


the  forms 


=  r„(-0. 114, 4.09  X  10"^)^  +  (0.305,0.111)^  +  (2.0,  0.6)A2A3,  (5.1a) 


^  =  r„(7.52  X  10-2,-7.74  X  lO'^)^  +  (0.265,0.117)^  +  (0.15, -0.15)AiA3, 
ati  dr\  adi 

(5.16) 

^  =  r„(-2.4  X  10-2,  -0.14)^  +  (0.14, 0.13)|^  +  (0.3, 0.4)AiA^,  (5.1c) 

for  the  resonant  triple  given  by  (3.3a)  and  (4.3).  As  already  mentioned,  the  effect 
of  having  four  critical  layers  within  the  flow  makes  the  calculation  of  the  interaction 
coefficients  very  sensitive  and  we  beUeve  that  the  values  of  these  quantities  quoted  in 
(5.1)  are  correct  to  within  about  ten  per  cent. 

Naturally,  the  evolution  of  the  wave  amplitudes  Ai,  A2  and  A3  is  very  dependent 
upon  the  initial  conditions  and  the  behaviours  of  these  quantities  may,  in  general,  only 
be  determined  by  a  numerical  approach.  We  do  not  study  this  facet  of  the  problem  here 
but,  instead,  draw  an  analogy  between  the  forms  of  (5.1)  and  those  derived  by  Smith  k. 
Stewart  (22)  in  their  investigation  of  resonant  triad  nonlinear  interaction  in  boundary 
layer  transition.  The  evolution  equations  derived  by  these  authors  were  a  simpler  form 
of  (5.1),  (effectively  with  dldr\  and  d/dOi  replaced  by  constants)  with  the  interaction 
coefficients  purely  imaginary  quantities.  Smith  k  Stewart  integrated  their  equations 
using  a  predictor-corrector  formulation  to  march  forward  in  time  and  the  resulting 
evolution  of  the  wave  amplitudes  took  a  very  complicated  form.  These  amplitudes 
generally  increased  with  time  although  this  trend  was  modified  by  the  superposition 
of  a  very  short  period,  large  oscillation  which  had  the  overall  effect  of  giving  the 
function  |Aj  (ti)|  a  very  spiky  appearance.  Similar  types  of  evolution  equations  have 
been  derived  by  Avis  (27),  who  was  concerned  with  resonant  triads  in  axisymmetric 
flow  on  a  cylindrical  geometry,  and  by  Bowles  (28),  who  investigated  triads  within 
transonic  flow.  These  authors  obtained  numerical  solutions  for  the  development  of  the 
mode  amplitudes,  and  both  found  behaviour  qualitative  to  that  described  in  (22). 

Very  little  is  known  about  particular  solutions  of  (5.1),  although  some  work  on 
related  equations  has  been  studied  by  Craik  (1,  29,  30).  No  general  analytic  solutions 
for  these  equations  have  been  found.  However,  an  important  property  of  the  system 
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(5.1)  is  under  what  conditions  these  equations  admit  singularities  in  the  solutions  Aj , 
(i  =  1,2,3),  either  in  a  finite  or  an  infinite  time.  Wang  (31)  has  derived  bounds  for 
the  solutions  of  three-wave  interaction  systems  with  arbitrary  coupling  coefficients. 
He  showed  that  if 

<^00  .  ,  dai  .  ,  da2  -  .  , 

-  iwoao  +  M0aia2)  "^  =  iwiai  + /iiaoa2,  =  ^0J2a2  +  M20oai,  (5.2) 

with  initial  data  ai(0)  =  OiO,  i  =  1,2,3;  then  explosive  instabilities  (instabilities  that 
become  unbounded  eventually,  possibly  after  infinite  time)  cannot  exist  if 


a  <  0  and  |ja(0)||  <  -\/|Q;|/y3, 

(5.3a) 

where 

a  =  min{  max{Aoo,  A12,  A21},  max{Aoi,  Aio,  A22}, 
max{Aoo,  Aoi,  A22},  max{Aoi,  A12,  A21}, 
max{Ao2,  Aio,  A21},  max{Ao2>  -^11,  A20}}, 

(5.35) 

/?  =  i  m^^{l/z<|},  and  Aj,  =  |/i,- 1  -  2Im  (w<) . 

(5.3c) 

Here  Hajj  = 

/  0  \  1/2 

ll(ao,ai,a2)ll  = 

If  we  allow  the  disturbances  proportional  to  (4.1)  to  be  neutrally  stable  in  space 
co-ordinates  ri  and  0i  we  can  substitute 


-  and  i0i,  (5.4) 

for  d(dri  and  djdQx  in  (5.1)  with  ai  and  real.  Then  using  the  results  (5.3)  we 
can  find  the  region  in  (ai,/?i)  space  for  which  explosive  instabilities  of  the  system 
(5.1)  are  impossible,  at  least  for  sufficiently  small  initial  disturbance  amplitudes.  This 
information  is  illustrated  in  Figure  (4)  where  we  indicate  the  minimum  size  of  the  initial 
amplitude  required  in  order  to  ensure  that  the  possibility  of  explosive  breakdown  is 
maintained.  Broadly  speaking,  we  can  conclude  from  figure  (4)  that  once  we  are 
located  in  that  part  of  (q;i,/9i)  space  at  which  a  finite  initial  amplitude  is  required  to 
enable  explosive  breakdown  then  by  decreasing  the  radial  wavenumber  ai  and/or  by 
increasing  the  azimuthal  wavenumber  /?i  there  needs  to  be  a  corresponding  rise  in  the 
initial  disturbance  amplitude  if  the  breakdown  mechanism  is  to  remain  possible. 
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It  should  be  noted  that  the  point  corresponding  to  the  choice  ai  =  =  0 

lies  in  the  zone  for  which  the  possibility  of  explosive  instabihty  cannot  be  discounted 
irrespective  of  the  initial  conditions.  Then  the  governing  equations  reduce  to  the  forms 

^  =  (2.0, 0.6)>l2A3,  ^  =  (0.15,  -0.15)^lA^,  ^  =  (0.3,  0.4) A1A2. 

^  ^  ^  (5.5) 

These  are,  of  course,  the  evolution  equations  obtained  by  allowing  the  disturbances 
to  evolve  temporally  whilst  suppressing  spatial  dependence.  It  is  a  straightforward  task 
to  write  61  =  A j,  62  —  .^2)  ^3  =  .^3  and  to  normalise  so  that  the  system  (5.5)  may  be 
rewritten  as 

dbi 

—  =  Sibjbl,  (i,i,fc  =  1,2,3  cyclically),  (5.6) 

where  the  interaction  coefficients  s,-  =  exp(i<^<).  Craik  (29,  30)  has  shown  that  ‘phase- 
locked’  solutions  of  systems  analogous  to  (5.6)  exist  where  0  =  +  ^2  +  ^3  is  constant 

and  Oi  =  phase{bi).  Then  0  must  be  a  root  of 

3 

tan(0<  —  0)  =  0. 

<=1 

In  this  case  the  existence  of  roots  of  this  equation  is  guaranteed  and  then,  for  ap¬ 
propriate  initial  conditions,  the  solutions  can  exhibit  a  finite  time  singularity.  It 
has  to  be  remembered  that  in  practice  once  the  disturbance  amplitudes  become  too 
large  however,  then  the  assumptions  concerning  the  smallness  of  the  amplitude  within 
definitions  (4.3)  cease  to  remain  valid  and  a  revised  asymptotic  analysis  is  required. 

Apart  from  the  results  relating  to  explosive  instabilities  described  above,  very  little 
is  known  about  solutions  of  system  (5.2)  for  general  parameter  values.  In  the  course 
of  an  investigation  into  interaction  processes  between  crossflow  vortices  in  practically 
important  cases  one  would  have  to  resort  to  extensive  numerical  calculations  to  provide 
a  full  description  of  the  characteristics  of  the  solution  of  the  evolution  equations.  This 
provides  scope  for  extending  the  elementary  ideas  presented  here  to  the  potentially 
more  important  interactions  in  experimentally  generated  boundary  layers. 

The  derivation  of  the  evolution  equations  relating  to  the  second  resonant  triad  as 
given  by  (3.36)  was  also  considered.  However,  for  this  triad,  the  mode  proportional  to 
Ez  has  a  critical  layer  very  close  to  the  wall.  In  practice,  it  was  found  that  the  numerical 
scheme  described  in  §4  was  too  crude  to  accurately  model  the  singular  behaviour  close 
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to  the  wall  to  provide  satisfactory  convergence.  This  discovery  points  to  the  need  to 
develop  improved  numerical  procedures  to  compute  the  interaction  coefficients. 

The  asymptotic  arguments  behind  the  derivation  of  the  amplitude  equations  for 
modes  with  spatial  and  temporal  dependences  as  in  (4.1),  (4.2)  relied  on  the  assump¬ 
tion  that  the  amplitudes  of  the  disturbances  were  0(S)  across  the  majority  of  the 
boundary  layer,  where  ^  1.  This  condition  was  introduced  primarily  to  ensure  that 

the  viscous  critical  layers  within  the  flow  assumed  their  simplest  possible  forms:  i.e. 
linear  critical  layers  encompassing  the  usual  ‘jtt’  phase  shift.  We  can  formally  relax 
this  requirement  by  making  the  following  observations.  Recall  that 


exp 


Oj  (r,  e)dr  dj3j  (i)  -  tflj 


(4.1a) 


and  suppose  that,  according  to  linear  theory,  the  mode  proportional  to  E,  is  neutrally 
stable  when 

Oj  =OiO}  +  eaij  -b . +  e*  d*,-  -t-  .  . . , 

Pi  =30j  +  ePij  + . -f-  Pk}  +  •  ■  • )  (5.7) 

Clj  =flOj  +  4" . +  +  •  •  •  • 

Suppose  we  allow  the  disturbance  to  evolve  on  slow  space  and  time  scales  given  by 
(4.2)  when  we  replace 

and  make  similar  replacements  for  6{d/dyi)  and  6(d/dti).  Then  the  analysis  of  the 
interaction  mechanism  will  proceed  in  a  very  similar  way  as  before,  so  long  as  the 
critical  layers  remain  linear  in  form.  Following  the  work  of  (8,  20),  we  know  that 
the  critical  layers  appearing  in  this  type  of  disturbance  structure  remain  linear  in 
character  until  the  disturbance  amplitude  reaches  tne  size  O(e^)  at  which  point  the 
critical  layer  configuration  is  modified  due  to  the  presence  of  (relatively)  large  mean 
flow  corrections  and  the  calculation  of  the  phase  shift  across  the  layer  becomes  a 
numerical  task.  However,  we  can  see  that  the  derivation  of  the  evolution  triad  of 
equations  (4.7)  may  be  trivially  modified  to  account  for  the  interaction  between  three 
modes  of  finite  amplitudes,  so  long  as  these  amplitudes  remain  less  than  O(e^). 

To  conclude,  we  believe  that  the  major  finding  of  our  work  here  is  that  triads 
of  linearly  neutrally  stable  crossflow  vortices  can  exist  in  three-dimensional  boundary 
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layers;  further  these  triads  are  such  that  \/e  have  a  resonant  set  of  modes  at  lead¬ 
ing  orders  in  wavenumber/  wavespeed  expansions  based  on  the  asymptotically  large 
Reynolds  number.  Further,  we  have  developed  a  fuUy  rational  description  of  the  in¬ 
teraction  processes  which  promises  to  reveal  the  degree  to  which  these  mechanisms 
can  alter  the  stability  characteristics  of  the  flow.  We  have  concentrated  on  a  special 
three-dimensional  flow  as  an  illustrative  example  but  the  results  obtained  here  suggest 
that  the  interaction  process  will  be  operative  in  more  practically  oriented  boundary 
layer  flows.  Our  work  can  be  regarded  as  a  first  stage  in  a  very  complicated  field: 
obvious  generalisations  could  be  tackled  to  explain  situations  in  which,  for  example, 
we  have  critical  layers  of  other  than  the  simple  linear  kind;  or  in  which  the  modes 
have  shght  detuning.  In  this  context  we  note  that  Avis  (27)  has  made  a  careful  study 
of  the  effects  of  detuning  on  resonant  triads  in  his  axisymmetric  flows.  The  analysis 
adopted  here  allows  nonparallel  effects  to  be  accounted  for  without  formal  difficulties- 
these  effects  can  play  significant  roles  in  practical  applications.  The  largest  drawback 
of  our  approach  is  that  the  computation  of  the  interaction  coefficients  appears  to  be 
not  straightforward,  especially  when  the  modes  have  critical  points  close  together  or 
when  such  critical  points  ate  near  the  bounding  surface.  We  feel  that  to  be  able  to  use 
our  approach  to  its  best  advantage  further  study  needs  to  be  devoted  to  developing 
more  powerful  and  sophisticated  techniques  for  resolution  of  the  numerical  problem 
as  opposed  to  the  fairly  rudimentary  ideas  employed  here.  Overall  however,  we  have 
here  the  first  asymptotic  description  by  which  the  importance  of  interaction  mecha¬ 
nisms  upon  a  resonant  triad  of  linear  neutrally  stable  crossflow  vortices  may  be  truly 
assessed. 
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Appendix. 


Here  we  briefly  describe  the  derivation  of  the  constants  ,  D^j }  appearing  in 
(4.6).  We  shall,  for  convenience,  determine  these  quantities  when  j  =  3  for  then  we 
may  easily  modify  the  results  for  j  —  2. 

With  /,  as  defined  in  (4.56),  we  suppose  that  close  to  the  critical  layer  ^3  the 
eigenfunction  Woo3  develops  according  to 

W^003  ~  ^3  +  O  ((^  -  6)  In  1^  -  61)  •  (Al) 

Further,  in  this  region  the  mode  proportional  to  E2  has  a  regular  Taylor  expansion 
and  so 

(U^002,V002,^002)  ^  (Qu02>Qv02,Qw02)  +  (^  -  Cs)  (Qu02>  ■5«02>  Qtu  02)  +•••»  (^2) 


where  the  quantities  Q|]o2}  constants.  Finally,  we  know  that  the  behaviours  of 
Uoo3  and  Vb03  can  be  obtained  by  suitable  manipulations  of  the  basic  equations  (2.6). 
Substituting  this  information  into  (4.56)  yie’ '  .:e  behaviour 


^  ‘o  Aia  X23 

(C  -  6)^  U  -  6)^  -  6) 


(A3) 


where 


^03  =  (aoi/003  -  Q!03/?0l)  , 

7tP03-O031 


Ai3  =  (Q:03/3o1  -  aoi/^os)  1(/?01  +  Pod)  Qv02  +  (ooi  +  0:03)  Qu02] 


and 


7t^03-Bo31 

ty|3-5032-B01lQv/02  _  jJ\Bo^ 
^031  7r/3o35o31 


.t 


(A4) 


(A5) 


(Q!01/?03  -  a03/?0l)  Qw02. 


i903Bo32Boi2Qw02  tQ«,027oi  1-^0311  ,  o  ^  \  t 

X23  =  - ^ - 'J  p-  - - (O(02P03  -  Q!03P02)  d- 


B, 


031 


7r/?03-Bo31 


— [(Q!0iQu02  +  /3oiQv02)  -Boil  +  (ao3/Joi  -  «0l/^03)  -  Qv02fnQj-J  ]  • 

■“031 

(A6) 


Here 


j  _  7rQ03Bo31fl'|3  _  ^^3 

"  703l-S03l|  |Bo3i| 


Then  we  can  show  that 


Wiai  ~ 


■Do3 
-  6) 


+  T)i3  In  -  .fa)  +  Di2  +  Dzz  U  -  -fa)  In  (^  -  ^3)  +  . . . , 


(^7) 


for  ^  >  ^3  with 


■D03  = 


^03 
2^010  ’ 


Diz  = 


2^011  ^"13 

■Bolo  -Boio  ’ 


and 


-D33  = 


^011 

Bqiq 


Di3  +  701^03  + 


X33 

■BolO 


The  quantities  Z)i3,  D33  appear  in  the  evolution  equa  ion  (4.7)  for  the  mode  pro¬ 
portional  to  £1.  Of  course,  we  need  to  obtain  quantities  equivalent  to  those  given 
above  in  order  to  traverse  the  critical  layer(s)  ^2>  but  these  can  be  deduced  merely  by 
interchanging  the  subscripts  ‘2’  and  ‘3’  in  the  relevant  formulae  above. 
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Figure  (1). 

Flow  structure  for  the  Rayleigh  modes  described  in  §2.  The  inviscid  layer  I  is  of 
thickness  O(f^)  (the  same  depth  as  the  boundary  layer  and  of  the  same  size  as  the 
horizontal  scale  for  the  structure).  In  addition,  we  have  a  viscous  wall  layer  (//)  of 
depth  and  a  critical  layer  (III)  of  thickness  O(e^)  at  those  positions  where 

the  ‘effective’  basic  velocity  profile  Vq  (see  (2.7))  vanishes. 


Figure  (2). 

b)  Illustration  that  for  0  <  p  <  /i,  there 
U,  =  CO  for  some  given  cq.  Hence  there  can 
flow  structure. 


Figure  (3). 

Solutions  of  the  eigenproblem  (2.7)  for  various  values  of  the  waveangle  parameter  M- 
a)  Radial  wavenumber  oq  as  a  function  of  the  azimuthal  wavenumber  /^o- 


Figure  (3). 

Solutions  of  the  eigenproblem  (2.7)  for  various  values  of  the  waveangle  parameter 
b)  Frequency  flo  of  mode  as  a  function  of  /Jo- 
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